Does stress mess with rodents' heads? Influence of habitat amount and genetic factors in mandible fluctuating asymmetry in South American water rats (Nectomys squamipes, Sigmodontinae) from Brazilian Atlantic rainforest remnants.

Loss of developmental stability can lead to deviations from bilateral symmetry (i.e. Fluctuating Asymmetry - FA), and is thought to be caused by environmental and genetic factors associated with habitat loss and stress. Therefore, levels of FA might be a valuable tool to monitor wild populations if FA serves as an indicator of exposure to stress due to impacts of habitat loss and fragmentation. In studies examining FA and habitat fragmentation, FA levels are often explained by loss of genetic variation, though few studies have addressed FA's use as indicator of environmental impact. Here, we investigated whether habitat loss, genetic variation, and/or inbreeding affect the developmental instability in Brazilian Atlantic forest populations of a Neotropical water rat (Nectomys squamipes). We sampled individuals from eight sites within Atlantic forest remnants with different amounts of available forest habitat and assessed FA levels with geometric morphometric techniques using adult mandibles. We used observed heterozygosity (Ho) and inbreeding coefficient (Fis), from seven microsatellite markers, as a proxy of genetic variation at individual and population levels. Populations were not significantly different for shape or size FA levels. Furthermore, interindividual variation in both shape and size FA levels and interpopulational differences in size FA levels were best explained by chance. However, habitat amount was negatively associated with both interpopulational variance and average shape FA levels. This association was stronger in populations living in areas with <28% of forest cover, which presented higher variance and higher average FA, suggesting that Nectomys squamipes might have a tolerance threshold to small availability of habitat. Our work is one of the first to use FA to address environmental stress caused by habitat loss in small mammal populations from a Neotropical biome. We suggest that shape FA might serve as a conservation tool to monitor human impact on natural animal populations.

Influence of female cuticular hydrocarbon (CHC) profile on male courtship behavior in two hybridizing field crickets Gryllus firmus and Gryllus pennsylvanicus.

BACKGROUND: The hybridizing field crickets, Gryllus firmus and Gryllus pennsylvanicus have several barriers that prevent gene flow between species. The behavioral pre-zygotic mating barrier, where males court conspecifics more intensely than heterospecifics, is important because by acting earlier in the life cycle it has the potential to prevent a larger fraction of hybridization. The mechanism behind such male mate preference is unknown. Here we investigate if the female cuticular hydrocarbon (CHC) profile could be the signal behind male courtship. RESULTS: While males of the two species display nearly identical CHC profiles, females have different, albeit overlapping profiles and some females (between 15 and 45%) of both species display a male-like profile distinct from profiles of typical females. We classified CHC females profile into three categories: G. firmus-like (F; including mainly G. firmus females), G. pennsylvanicus-like (P; including mainly G. pennsylvanicus females), and male-like (ML; including females of both species). Gryllus firmus males courted ML and F females more often and faster than they courted P females (p < 0.05). Gryllus pennsylvanicus males were slower to court than G. firmus males, but courted ML females more often (p < 0.05) than their own conspecific P females (no difference between P and F). Both males courted heterospecific ML females more often than other heterospecific females (p < 0.05, significant only for G. firmus males). CONCLUSIONS: Our results suggest that male mate preference is at least partially informed by female CHC profile and that ML females elicit high courtship behavior in both species. Since ML females exist in both species and are preferred over other heterospecific females, it is likely that this female type is responsible for most hybrid offspring production.

Mapping reduced introgression loci to the X chromosome of the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus.

The genomic architecture of barriers to gene exchange during the speciation process is poorly understood. The genomic islands model suggests that loci associated with barriers to gene exchange prevent introgression of nearby genomic regions via linkage disequilibrium. But few analyses of the actual genomic location of non-introgressing loci in closely related species exist. In a previous study Maroja et al. showed that in the hybridizing field crickets, Gryllus firmus and G. pennsylvanicus, 50 non-introgressing loci are localized on two autosomal regions and the X chromosome, but they were not able to map the loci along the X chromosome because they used a male informative cross. Here, we localize the introgressing and non-introgressing loci on the X chromosome, and reveal that all X-linked non-introgressing loci are restricted to a 50-cM region with 10 of these loci mapped to a single location. We discuss the implications of this finding to speciation.

Rapid sexual and genomic isolation in sympatric Drosophila without reproductive character displacement.

The rapid evolution of sexual isolation in sympatry has long been associated with reinforcement (i.e., selection to avoid maladaptive hybridization). However, there are many species pairs in sympatry that have evolved rapid sexual isolation without known costs to hybridization. A major unresolved question is what evolutionary processes are involved in driving rapid speciation in such cases. Here, we focus on one such system; the Drosophila athabasca species complex, which is composed of three partially sympatric and interfertile semispecies: WN, EA, and EB. To study speciation in this species complex, we assayed sexual and genomic isolation within and between these semispecies in both sympatric and allopatric populations. First, we found no evidence of reproductive character displacement (RCD) in sympatric zones compared to distant allopatry. Instead, semispecies were virtually completely sexually isolated from each other across their entire ranges. Moreover, using spatial approaches and coalescent demographic simulations, we detected either zero or only weak heterospecific gene flow in sympatry. In contrast, within each semispecies we found only random mating and little population genetic structure, except between highly geographically distant populations. Finally, we determined that speciation in this system is at least an order of magnitude older than previously assumed, with WN diverging first, around 200K years ago, and EA and EB diverging 100K years ago. In total, these results suggest that these semispecies should be given full species status and we adopt new nomenclature: WN-D. athabasca, EA-D. mahican, and EB-D. lenape. While the lack of RCD in sympatry and interfertility do not support reinforcement, we discuss what additional evidence is needed to further decipher the mechanisms that caused rapid speciation in this species complex.

The wavy Mutation Maps to the Inositol 1,4,5-Trisphosphate 3-Kinase 2 (IP3K2) Gene of Drosophila and Interacts with IP3R to Affect Wing Development.

Inositol 1,4,5-trisphosphate (IP3) regulates a host of biological processes from egg activation to cell death. When IP3-specific receptors (IP3Rs) bind to IP3, they release calcium from the ER into the cytoplasm, triggering a variety of cell type- and developmental stage-specific responses. Alternatively, inositol polyphosphate kinases can phosphorylate IP3; this limits IP3R activation by reducing IP3 levels, and also generates new signaling molecules altogether. These divergent pathways draw from the same IP3 pool yet cause very different cellular responses. Therefore, controlling the relative rates of IP3R activation vs. phosphorylation of IP3 is essential for proper cell functioning. Establishing a model system that sensitively reports the net output of IP3 signaling is crucial for identifying the controlling genes. Here we report that mutant alleles of wavy (wy), a classic locus of the fruit fly Drosophila melanogaster, map to IP3 3-kinase 2 (IP3K2), a member of the inositol polyphosphate kinase gene family. Mutations in wy disrupt wing structure in a highly specific pattern. RNAi experiments using GAL4 and GAL80(ts) indicated that IP3K2 function is required in the wing discs of early pupae for normal wing development. Gradations in the severity of the wy phenotype provide high-resolution readouts of IP3K2 function and of overall IP3 signaling, giving this system strong potential as a model for further study of the IP3 signaling network. In proof of concept, a dominant modifier screen revealed that mutations in IP3R strongly suppress the wy phenotype, suggesting that the wy phenotype results from reduced IP4 levels, and/or excessive IP3R signaling.

Genes with Restricted Introgression in a Field Cricket (Gryllus firmus/Gryllus pennsylvanicus) Hybrid Zone Are Concentrated on the X Chromosome and a Single Autosome.

Characterizing the extent of genomic differentiation between recently diverged lineages provides an important context for understanding the early stages of speciation. When such lineages form discrete hybrid zones, patterns of differential introgression allow direct estimates of which genome regions are likely involved in speciation and local adaptation. Here we use a backcross experimental design to construct a genetic linkage map for the field crickets Gryllus firmus and Gryllus pennsylvanicus, which interact in a well-characterized hybrid zone in eastern North America. We demonstrate that loci with major allele frequency differences between allopatric populations are not randomly distributed across the genome. Instead, most are either X-linked or map to a few small autosomal regions. Furthermore, the subset of those highly differentiated markers that exhibit restricted introgression across the cricket hybrid zone are also concentrated on the X chromosome (39 of 50 loci) and in a single 7-cM region of one autosome. Although the accumulation on the sex chromosome of genes responsible for postzygotic barriers is a well-known phenomenon, less attention has been given to the genomic distribution of genes responsible for prezygotic barriers. We discuss the implications of our results for speciation, both in the context of the role of sex chromosomes and also with respect to the likely causes of heterogeneous genomic divergence. Although we do not yet have direct evidence for the accumulation of ecological, behavioral, or fertilization prezygotic barrier genes on the X chromosome, faster-X evolution could make these barriers more likely to be X-linked.

Development and characterization of 10 microsatellite markers in Sagina nodosa (Caryophyllaceae).

PREMISE OF THE STUDY: We developed 10 novel microsatellite loci for Sagina nodosa, a diploid perennial arctic-alpine herb. To our knowledge, these are the first microsatellite loci for a Sagina species. • METHODS AND RESULTS: We performed a low-coverage 454 next-generation sequencing of enriched genomic fragments derived from one individual to generate a massive library of contigs containing potential polymorphic microsatellites. We present data for 10 novel polymorphic microsatellite loci containing di-, tri-, tetra-, and hexanucleotide repeats with two to nine alleles per locus assessed in 29 individuals. • CONCLUSIONS: These polymorphic microsatellite loci in S. nodosa will provide insights on the population structure and life history of S. nodosa in Isle Royale and other North American populations.

Barriers to gene exchange in hybridizing field crickets: the role of male courtship effort and cuticular hydrocarbons.

BACKGROUND: Pre-zygotic barriers often involve some form of sexual selection, usually interpreted as female choice, as females are typically the choosier sex. However, males typically show some mate preferences, which are increasingly reported. Here we document previously uncharacterized male courtship behavior (effort and song) and cuticular hydrocarbon (CHC) profiles in the hybridizing crickets Gryllus firmus and G. pennsylvanicus. These two species exhibit multiple barriers to gene exchange that act throughout their life history, including a behavioral barrier that results in increased time to mate in heterospecific pairs. RESULTS: We demonstrated that male mate choice (as courtship effort allocation) plays a more important role in the prezygotic behavioral barrier than previously recognized. In gryllids females ultimately decide whether or not to mate, yet we found males were selective by regulating courtship effort intensity toward the preferred (conspecific) females. Females were also selective by mating with more intensely courting males, which happened to be conspecifics. We report no differences in courtship song between the two species and suggest that the mechanism that allows males to act differentially towards conspecific and heterospecific females is the cuticular hydrocarbon (CHC) composition. CHC profiles differed between males and females of both species, and there were clear differences in CHC composition between female G. firmus and G. pennsylvanicus but not between the males of each species. CONCLUSION: Although many barriers to gene exchange are known in this system, the mechanism behind the mate recognition leading to reduced heterospecific mating remains unknown. The CHC profiles might be the phenotypic cue that allow males to identify conspecifics and thus to adjust their courtship intensity accordingly, leading to differential mating between species.

High-throughput microsatellite marker development for the distylous herb Primula mistassinica (Primulaceae).

PREMISE OF THE STUDY: Twelve microsatellite markers were developed for Primula mistassinica, a distylous, diploid arctic-alpine plant. The markers will be used to investigate the landscape genetics of a disjunct population on Isle Royale, Michigan, and the phylogeographic patterns of the species. • METHODS AND RESULTS: We used Roche/454 high-throughput technology to sequence microsatellite-enriched regions in the P. mistassinica genome. We developed 12 polymorphic microsatellite primer sets. These loci contained di-, tri-, and tetranucleotide repeats with two to nine alleles per locus when assessed in 23 individuals. • CONCLUSIONS: Understanding the historical movements of P. mistassinica will provide insight to the survival prospects of current Arctic plant populations, which face the pressures of global, anthropogenic climate change.

Partial complementarity of the mimetic yellow bar phenotype in Heliconius butterflies.

Heliconius butterflies are an excellent system for understanding the genetic basis of phenotypic change. Here we document surprising diversity in the genetic control of a common phenotype. Two disjunct H. erato populations have each recruited the Cr and/or Sd loci that control similar yellow hindwing patterns, but the alleles involved partially complement one another indicating either multiple origins for the patterning alleles or developmental drift in genetic control of similar patterns. We show that in these H. erato populations cr and sd are epistatically interacting and that the parental origin of alleles can explain phenotypes of backcross individuals. In contrast, mimetic H. melpomene populations with identical phenotypes (H. m. rosina and H. m. amaryllis) do not show genetic complementation (F(1)s and F(2)s are phenotypically identical to parentals). Finally, we report hybrid female inviability in H. m. melpomene × H. m. rosina crosses (previously only female infertility had been reported) and presence of standing genetic variation for alternative color alleles at the Yb locus in true breeding H. melpomene melpomene populations (expressed when in a different genomic background) that could be an important source of variation for the evolution of novel phenotypes or a result of developmental drift. Although recent work has emphasized the simple genetic control of wing pattern in Heliconius, we show there is underlying complexity in the allelic variation and epistatic interactions between major patterning loci.

Characterisation and expression of microRNAs in developing wings of the neotropical butterfly Heliconius melpomene.

BACKGROUND: Heliconius butterflies are an excellent system for studies of adaptive convergent and divergent phenotypic traits. Wing colour patterns are used as signals to both predators and potential mates and are inherited in a Mendelian manner. The underlying genetic mechanisms of pattern formation have been studied for many years and shed light on broad issues, such as the repeatability of evolution. In Heliconius melpomene, the yellow hindwing bar is controlled by the HmYb locus. MicroRNAs (miRNAs) are important post-transcriptional regulators of gene expression that have key roles in many biological processes, including development. miRNAs could act as regulators of genes involved in wing development, patterning and pigmentation. For this reason we characterised miRNAs in developing butterfly wings and examined differences in their expression between colour pattern races. RESULTS: We sequenced small RNA libraries from two colour pattern races and detected 142 Heliconius miRNAs with homology to others found in miRBase. Several highly abundant miRNAs were differentially represented in the libraries between colour pattern races. These candidates were tested further using Northern blots, showing that differences in expression were primarily due to developmental stage rather than colour pattern. Assembly of sequenced reads to the HmYb region identified hme-miR-193 and hme-miR-2788; located 2380 bp apart in an intergenic region. These two miRNAs are expressed in wings and show an upregulation between 24 and 72 hours post-pupation, indicating a potential role in butterfly wing development. A search for miRNAs in all available H. melpomene BAC sequences (~2.5 Mb) did not reveal any other miRNAs and no novel miRNAs were predicted. CONCLUSIONS: Here we describe the first butterfly miRNAs and characterise their expression in developing wings. Some show differences in expression across developing pupal stages and may have important functions in butterfly wing development. Two miRNAs were located in the HmYb region and were expressed in developing pupal wings. Future work will examine the expression of these miRNAs in different colour pattern races and identify miRNA targets among wing patterning genes.

Genomic hotspots for adaptation: the population genetics of Müllerian mimicry in the Heliconius melpomene clade.

Wing patterning in Heliconius butterflies is a longstanding example of both Müllerian mimicry and phenotypic radiation under strong natural selection. The loci controlling such patterns are "hotspots" for adaptive evolution with great allelic diversity across different species in the genus. We characterise nucleotide variation, genotype-by-phenotype associations, linkage disequilibrium, and candidate gene expression at two loci and across multiple hybrid zones in Heliconius melpomene and relatives. Alleles at HmB control the presence or absence of the red forewing band, while alleles at HmYb control the yellow hindwing bar. Across HmYb two regions, separated by approximately 100 kb, show significant genotype-by-phenotype associations that are replicated across independent hybrid zones. In contrast, at HmB a single peak of association indicates the likely position of functional sites at three genes, encoding a kinesin, a G-protein coupled receptor, and an mRNA splicing factor. At both HmYb and HmB there is evidence for enhanced linkage disequilibrium (LD) between associated sites separated by up to 14 kb, suggesting that multiple sites are under selection. However, there was no evidence for reduced variation or deviations from neutrality that might indicate a recent selective sweep, consistent with these alleles being relatively old. Of the three genes showing an association with the HmB locus, the kinesin shows differences in wing disc expression between races that are replicated in the co-mimic, Heliconius erato, providing striking evidence for parallel changes in gene expression between Müllerian co-mimics. Wing patterning loci in Heliconius melpomene therefore show a haplotype structure maintained by selection, but no evidence for a recent selective sweep. The complex genetic pattern contrasts with the simple genetic basis of many adaptive traits studied previously, but may provide a better model for most adaptation in natural populations that has arisen over millions rather than tens of years.

Genealogical discordance and patterns of introgression and selection across a cricket hybrid zone.

In recently diverged species, ancestral polymorphism and introgression can cause incongruence between gene and species trees. In the face of hybridization, few genomic regions may exhibit reciprocal monophyly, and these regions, usually evolving rapidly under selection, may be important for the maintenance of species boundaries. In animals with internal fertilization, genes encoding seminal protein are candidate barrier genes. Recently diverged hybridizing species such as the field crickets Gryllus firmus and G. pennsylvanicus, offer excellent opportunities to investigate the origins of barriers to gene exchange. These recently diverged species form a well-characterized hybrid zone, and share ancestral polymorphisms across the genome. We analyzed DNA sequence divergence for seminal protein loci, housekeeping loci, and mtDNA, using a combination of analytical approaches and extensive sampling across both species and the hybrid zone. We report discordant genealogical patterns and differential introgression rates across the genome. The most dramatic outliers, showing near-zero introgression and more structured species trees, are also the only two seminal protein loci under selection. These are candidate barrier genes with possible reproductive functions. We also use genealogical data to examine the demographic history of the field crickets and the current structure of the hybrid zone.

Searching for candidate speciation genes using a proteomic approach: seminal proteins in field crickets.

In many animals, male seminal proteins influence gamete interactions and fertilization ability and are probably involved in barriers to gene flow between diverging lineages. Here we use a proteomic approach to identify seminal proteins that are transferred to females during copulation and that may be involved in fertilization barriers between two hybridizing field crickets (Gryllus firmus and Gryllus pennsylvanicus). Analyses of patterns of divergence suggest that much of the field cricket genome has remained undifferentiated following the evolution of reproductive isolation. By contrast, seminal protein genes are highly differentiated. Tests of selection reveal that positive selection is likely to be responsible for patterns of differentiation. Together, our observations suggest that some of the loci encoding seminal proteins may indeed play a role in fertilization barriers in field crickets.

Phylogeography of spruce beetles (Dendroctonus rufipennis Kirby) (Curculionidae: Scolytinae) in North America.

Tree-feeding insects that are widespread in north temperate regions are excellent models for studying how past glaciations have impacted differentiation and speciation. We used mitochondrial DNA (mtDNA) sequences and allele frequencies at nine microsatellite loci to examine genetic population structure across the current range of the spruce beetle (Dendroctonus rufipennis), an economically important insect in North America. Two major haplotype groups occur across northern North America, from Newfoundland to Alaska, on white spruce (Picea glauca), and a third distinctive haplotype group occurs throughout the Rocky Mountains on Engelmann spruce (Picea engelmannii). The two mtDNA lineages found in northern populations are 3-4% divergent from each other and from the lineages found in the Rocky Mountains. Analyses of microsatellite data also suggest the existence of major population groupings associated with different geographical regions. In the Pacific Northwest, concordant contact zones for genetically distinct populations of spruce beetles and their principal hosts appear to reflect recent secondary contact. Although we could detect no evidence of historical mtDNA gene flow between allopatric population groups, patterns of variation in the Pacific Northwest suggest recent hybridization and introgression. Together with the pollen record for spruce, they also suggest that beetles have spread from at least three glacial refugia. A minimum estimate of divergence time between the Rocky Mountain and northern populations was 1.7 Myr (million years), presumably reflecting the combined effects of isolation during multiple glacial cycles.

Identification and comparative analysis of accessory gland proteins in Orthoptera.

Accessory reproductive gland proteins (Acps) in Drosophila evolve quickly and appear to play an important role in ensuring the fertilization success of males. Moreover, Acps are thought to be involved in establishing barriers to fertilization between closely related species. While accessory glands are known to occur in the males of many insect groups, the proteins that are passed on to females by males during mating have not been well characterized outside of Drosophila. To gain a better understanding of these proteins, we characterized ESTs from the accessory glands of two cricket species, Allonemobius fasciatus and Gryllus firmus. Using an expressed sequence tag (EST) approach, followed by bioinformatic and evolutionary analyses, we found that many proteins are secreted and, therefore, available for transfer to the female during mating. Further, we found that most ESTs are novel, showing little sequence similarity between taxa. Evolutionary analyses suggest that cricket proteins are subject to diversifying selection and indicate that Allonemobius is much less polymorphic than Gryllus. Despite rapid nucleotide sequence divergence, there appears to be functional conservation of protein classes among Drosophila and cricket taxa.

Molecular evolution of seminal proteins in field crickets.

In sexually reproducing organisms, male ejaculates are complex traits that are potentially subject to many different selection pressures. Recent experimental evidence supports the hypothesis that postmating sexual selection, and particularly sexual conflict, may play a key role in the evolution of the proteinaceous components of ejaculates. However, this evidence is based almost entirely on the study of Drosophila, a species with a mating system characterized by a high cost of mating for females. In this paper, we broaden our understanding of the role of selection on the evolution of seminal proteins by characterizing these proteins in field crickets, a group of insects in which females appear to benefit from mating multiply. We have used an experimental protocol that can be applied to other organisms for which complete genome sequences are not yet available. By combining an evolutionary expressed sequence tag screen of the male accessory gland in 2 focal species (Gryllus firmus and Gryllus pennsylvanicus) with a bioinformatics approach, we have been able to identify as many as 30 seminal proteins. Evolutionary analyses among 5 species of the genus Gryllus suggest that seminal protein genes evolve more rapidly than genes encoding proteins that are not involved with reproduction. The rates of synonymous substitution (dS) are similar in genes encoding seminal proteins and genes encoding "housekeeping" proteins. For the same comparison, the rate of fixation of nonsynonymous substitutions (dN) is 3 times higher in genes encoding seminal proteins, suggesting that the divergence of seminal proteins in field crickets has been accelerated by positive Darwinian selection. In spite of the contrasting characteristics of the Drosophila and Gryllus mating systems, the mean selection parameter omega and the proportion of loci estimated to be affected by positive selection are very similar.

Population structure and genetic variability of mainland and insular populations of the Neotropical water rat, Nectomys squamipes (Rodentia, Sigmodontinae)

Seven microsatellite loci were used to investigate the genetic variability and structure of six mainland and two island populations of the Neotropical water rat Nectomys squamipes, a South American semi-aquatic rodent species with a wide distribution. High levels of variability were found within mainland populations while island populations were less variable but the more differentiated in respect to allele number and frequency. The time of biological divergence between mainland and island populations coincided with geological data. A significant geographic structure was found in mainland populations (q = 0.099; r = 0.086) although the degree of differentiation was relatively low in respect to the distance between surveyed localities (24 to 740 km). Genetic and geographic distances were not positively correlated as previously found with random amplified polymorphic DNA (RAPD) markers. Significant but low genetic differentiation in the mainland and lack of isolation by distance can be explained by large population size and/or recent population expansion. Additionally, the agreement between the age of geologic events (sea level fluctuations) and divergence times for insular populations points to a good reference for molecular clock calibration to associate recent environmental changes and the distribution pattern of small mammals in the Brazilian Atlantic Forest